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ABSTRACT 

Context. The photospheric abundances of evolved solar-type stars of different metallicities serve as probes into stellar evolution the- 
ory. 

Aims. Stellar photospheres of post-asymptotic giant branch (post-AGB) stars bear witness to the internal chemical enrichment pro- 
cesses, integrated over their entire stellar evolution. Here we study post-AGB stars in the Large Magellanic Cloud (LMC). With their 
known distances, these rare objects are ideal tracers of AGB nucleosynthesis and dredge-up phenomena. 

Methods. We used the UVES spectrograph mounted on the Very Large Telescope (VLT) at the European Southern Observatory 
(ESO), to obtain high-resolution spectra with high signal-to-noise (S/N) of a sample of four post-AGB stars. The objects display a 
spectral energy distribution (SED) that indicates the presence of circumstellar dust. We perform a detailed abundance analysis on the 
basis of these spectra. 

Results. All objects are C-rich, and strongly enhanced in s-process elements. We deduced abundances of heavy s-process elements 
for all stars in the sample, and even found an indication of the presence of Hg in the spectrum of one object. The metallicity of all 
stars except J053253. 5 1-695915.1 is considerably lower than the average value that is observed for the LMC. The derived luminosities 
show that we witness the late evolution of low-mass stars with initial masses close to 1 M Q . An exception is J053253. 51-695915.1 
and we argue that this object is likely a binary. 

Conclusions. We confirmed the correlation between the efficiency of the third-dredge up and the neutron exposure that is detected 
in Galactic post-AGB stars. The non-existence of a correlation between metallicity and neutron irradiation is also confirmed and 
expanded to smaller metallicities. We confirm the status of 21 /im stars as post-Carbon stars. Current theoretical AGB models over- 
estimate the observed C/O ratios and fail to reproduce the variety of s-process abundance patterns that is observed in otherwise very 
similar objects. Similar results have recently been found for a post-AGB star in the Small Magellanic Cloud (SMC). 

Key words. Stars: AGB and post-AGB - Stars: evolution - (Galaxies:) Magellanic Clouds - (Stars:) circumstellar matter 



1. Introduction 

Post-AGB stars are stars of low and intermediate initial mass 
(M ^ 7) M Q ) that approach the end of their stellar evolution. 
The post-AGB evolutionary phase starts when the strong dusty 
mass loss, that is thought to characterise the end of the AGB, 
has ceased. The central star will evolve rapidly to the blue of the 
HR diagram. Here we concentrate on the photospheric chemi- 
cal abundances and want to use the photosphere of post-AGB 
stars as ideal tracers to study AGB nucleosynthetic and dredge- 
up processes. 

Among the known Galactic post-AGB stars, the variety in 
observed chemical patterns in the photospheres is large (e.g. 
Van Winc kel I2003T and references therein). Some objects are 
strongly enhanced in neutron capture (s-process) elements, while 
others are not enhanced at all and might even display s-process 
deficiencies. This bifurcation is rather strict, as only a limited 
number of post-AGB objects are known to be m ildly enriched 
in s-process elements (e.g. [Sumangala Rao et al]|2011l and ref- 
erences therein). This large chemical variety and the rather small 
sample of known Galactic post-AGB stars, makes it hard to un- 
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derstand how these results can be interpret in the framework of 
chemical AGB evolution models. 

One of the problems with the current sample of known 
Galactic post-AGB stars is that their distances, and hence lumi- 
nosities and initial masses, are badly determined as practically 
none of these objects has a reliable parallax measurement. We 
want to overcome this problem by studying post-AGB objects 
in other galaxies, and therefore focus in this contribution on the 
LMC. At a distance of about 50 kpc (Stor m et al.ll201 ll) and a t 
a favourable aspect angle of 35° ( van der Marel & Cioni 2001 ), 
the LMC is far enough away that all objects can be assumed to 
reside at the same distance, and still close enough to allow the 
observation of individual objects in detail. 

Observations of post-AGB stars can help us interpret the 
AGB phase of which they are the progeny. Studying optically 
bright post-AGB stars holds several advantages over studying 
the AGB stars themselves. First, a simultaneous study of the 
photosphere and the circumstellar environment is possible when 
the circumstellar envelope becomes optically thin by expan- 
sion. Photospheric spectra of AGB stars are veiled by molecu- 
lar lines which impedes spectral analysis of trace elements (e.g. 
lAbia et all 12008 ). Spectra of post-AGB stars are dominated by 
atomic transitions rather than molecular transitions and this al- 
lows abundance determinations of a wide range of elements. 
Furthermore, the photospheres of post-AGB stars are convenient 
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Tabl e 1. Names and spectral types of the programme stars 
(Ivan Aarle etal .11201 ll) . 



Table 2. Observational log. 



Object Name 


Spectral Type" 


(IRAC) 




J050632. 10-7 14229.8 


A3Iab/F0Ibp(e) 


J052043. 86-692341.0 


F5Ib(e) 


J053250.69-7 13925.8 


F6Ia 


J053253.51-695915.1 


G5 



Notes. (a) The '/' indicates that more than one low resolution, optical 
spectrum was obtained. 



to study as the spectral types tend to cluster in a range where 
model atmosphere analysis in the optical is the most reliable. 

We initiated a large project to study quantitatively our 
newly identified , large sample of post-AGB stars in the LMC 
(Ivan Aarle et all 1201 ll) and Small Magellanic Cloud (SMC) 
(Kamath etal. 2011). In this paper, we will present the abun- 
dance analysis of four LMC post-AGB stars that were selected 
from the 70 spectroscopically confirmed post-AGB stars in the 
catalogue of optical ly bright post-AGB stars of the LMC (see 
Ivan Aarle et al.ll201 ll) . We introduce the objects in Sect. [2] The 
high-resolution UVES spectra and their reduction is discussed in 
Sect. [3] We performed a detailed abundance analysis of the ob- 
tained spectra which is presented in Sect. [4] The results of this 
analysis are discussed in the next sections: In Sect. |6]we use the 
acquired atmosphere models to constrain the luminosities and 
initial masses of our program stars, the neutron exposure of the 
present sample is compared to that of some Galactic post-AGB 
objects in Sect. [7J Sect. [8] deals with the possibility that one of 
the objects in the sample is a binary. A first comparison with 
chemical AGB models is presented in Sect. [9] We end with some 
conclusions and future prospects in Sect. [10] 



2. The sample stars 

Basic information on four objects is gathered in Table Q] This 
table lists the full name of each object, which we will abbreviate 
in the rest of this paper to increase the readability. 

Three objects have a clearly double peaked SED (Fig. [7): in 
which the peak at longer wavelengths is indicative of a freely 
expanding, detached shell (van Aarle et al. 201 1). In the SED of 
J053253, the fourth object, only a small, relatively hot infrared 
excess is visible. 

All objects in our sample display photometric variability. 
J052043 and J053250 are semi-regular variables (SRVs) with 
an unusual warm spectral type, while J050632 and J053253 are 
population II Cepheids with a very low amplitude and pe riods 
of 49.1 and 91.04 days respectively ( Ivan Aarle et al.ll201 ll) . All 
objects have a peak-to-peak photometric amplitude of less than 
0.1 mag. 

Only one ob ject in our sample has been studied before. 
Volk et al. (201 1) obtained Spitzer Space Telescope mid-infrared 
spectroscopy for J052043 and detected a strong 21 jj.m fea- 
ture. The carrier of this dust feature is still unknown but i t 
is only detected in post-AGB stars (e.g. iHrivnaket al. 2009). 
For the Galactic sample of post-AGB stars, it is known that 
the presence of this feature accompanies a central star which 
is enriched in carbon and shows a strong enhancement in s- 
process elements ( Van Winckel & Revniers 2000; Reddy et al. 
2002; Reyniers 2002; Reyniers et al. 2004). In this paper, we will 
explore if this remains true for this object in the LMC. 



Object 


Date 


Exp. Time 




S/N 








(s) 


Blue 


Rediow 


Red up 


J050632 


4-5/11/2008 


2 x 3005 


50 


105 


110 


J052043 


21-22/11/2008 


1 x 3005 


10 


50 


60 




4-5/12/2008" 


2 x 3005 


15 










3 x 3005 




80 


100 


J053250 


17-18/11/2008 


3 x 3005 


15 


90 


120 


J053253 


19-20/10/2008 


4 x 3005 


10 


50 


70 



Notes. We list the names, dates, exposure times and signal-to-noise per 

pixel of all UVES spectra taken for our sample stars. 'Blue' corresponds 

to the spectral range from 3280 - 4519 A, 'Red low ' to 4777 - 5757 A, 

and 'Red up ' to 5834 - 6810 A. 

<o) The third blue spectrum had to be removed due to an instrumental 

error. 



3. Observations 

We obtained high-resolution, high-signal-to-noise UVES spec- 
tra (Dekker et al. 2000) for all stars in our sample. The spectra 
were taken in service mode in the fall of 2008. The spectrograph 
splits the light beam from the telescope into two spectral arms 
which are operated simultaneously. This results in a spectral cov- 
erage from 3280 to 4519 A in the Blue arm, and from 4777 to 
5757 A, and 5834 to 6810 A in the Red arm. We used a slit 
width of 0.8" which translates into 6 A ~ A/52 000 in the blue 
and 8 A ~ /1/48 000 in the red spectral regiorQ. Some details of 
the observations can be found in Table [2] 

The data reduction was performed with the UVES pipeline 
(version 4.3.0), which included bias correction, flat-fielding, and 
background, sky, and cosmic hit correction. Even though the im- 
age sheer was used to obtain our spectra, we used optimal in- 
stead of average extraction to convert the frames from pixel-pixel 
to pixel-order space. This produced spectra of similar quality but 
took better care of the remaining cosmic hits. 

Before merging the different spectra, the barycentric correc- 
tion of the spectra, as well as an intrinsic variability check on 
the different exposures of each object were performed. All spec- 
tra are consistent for all objects, except for J052043. Its spectra 
obtained in December show a small difference in spectral lines 
with respect to the ones taken in November. We decided not to 
merge these spectra, but to treat them separately. We will in- 
dicate the spectrum taken in November with a, and that from 
December with b. We removed some additional cosmic hits by 
hand, and merged all consistent spectra, thus retaining two spec- 
tra for J052043 and one for all other objects. All spectra were 
normalised by dividing the spectrum by a spline function de- 
fined through interactively identified continuum points. The re- 
ciprocal of the standard deviation computed on a few normalised 
continuum windows was used to express the S/N of our spectra 
(Tabled and some sample spectra are shown in Fig.Q] 

The high resolution and the large wavelength coverage of our 
UVES spectra permit us to determine the radial velocities using 
the position of a large number of lines. The line positions were 
measured using a Gaussian fit through the measured profiles. 
The identification of th e lines and their re st wavelength were 
taken from the VALD dKupkaet alj|1999l) database. The line 
identification itself was performed by comparing the expected 
line strength of all input lines using the model atmosphere and 
overabundances (Sect. 4) of our chemical analyses. We took au- 
tomatically the nearest strongest line. The internal consistancy 



http : //www . eso . org/sci/f acilities/paranal/instruments/uves 
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Fig. 1. Example of the quality of the spect ra of the stars in o ur sample. A complete line identification of this spectral interval has 
been performed using the VALD database (Kup ka et al.ll 19991) . Elements with atomic number Z smaller than 30 are indicated with a 
full line, light s-process elements with a dashed line, and heavy s-process elements with a dotted line. As the lower spectrum differs 
significantly from the upper four, other features are indicated on it. All objects are clearly enhanced in s-process elements, although 
this is less obvious in the spectrum of J050632 due to its higher temperature. 



Table 3. Heliocentric radial velocities. 



Object 


V,- 

(km/s) 


Nii„ cs 


J050632 


218.2 ±0.8 


364 


J052043_a 


243.3 ± 1.1 


346 


J052043_b 


244.6 ± 1.0 


438 


J053250 


283.3 ± 0.7 


278 


J053253 


237.5 ± 0.8 


454 



Notes. We list the names, radial velocities and the number of lines used 
to determine these velocities. The errors on the radial velocity are the 
standard deviations of the individual results of all lines. 



on a large number of lines shows that this method is reliable. 
The results are listed in Table[3] The radial velocities we find are 
consistent with the average radial velocity of the LMC, which is 
about 262 km/s ( van der Marel et al. 2002). 



4. Spectral Analyses 

We use the list of lines that are useful for the chemical analy- 
sis of A, F, and G stars that has been gathered at the Instituut 
voor Sterrenk unde and is regularly u pdated. The initial list is 
described in lVan Winckel & Revniersl (l2000h . and contains lines 
for a range of elements from He (Z = 2) up to Eu (Z = 63). 
With this list, we ensure that the different chemical analyses 
are performed homogeneously and we restrict ourselves to lines 
with reliable oscillator strengths. We determined the equivalent 
widths of the lines by fitting Gaussian profiles to the observa- 
tions to take care of blending in their wings. 



To determine the atmosphere parameters and obtain the 
abundances of trace elements, we u sed the ATLAS 9 model at- 
mospheres dCastelli & Kurucz 2004) in combination with the 
latest version (August 2010) of the LTE abundance calcula- 
tion routine MOOG dSnedenlll973l) . The model gridsteps are 
Ar eff = +/ - 250K, Alogg = +/ - 0.5 and A([Fe/H]) = +/- 
0.5 



4.1. Atmospheric parameters 

We used several spectroscopic tracers to quantify the effective 
temperature of our low gravity objects. A useful upper limit for 
the effective temperature is provided by the absence of He lines 
in the spectra. At the S/N of our spectra, this means that the lines 
must be weaker than ~ 5 mA. The absence of He lines, is espe- 
cially constraining for the hottest object in our sample J050632. 
In Table |U we list for several models the expected equivalent 
width of the strongest He features when assuming a solar he- 
lium abundance. The solar He abundance is a lower limit for the 
actual abundance as an increase of the helium abundance in the 
envelope of our stars of up to 20 % can be exp ected during the 
AGB evolution (e.g JCristallo et al.ll2009l.l201 ll) . 

An excitation analyses of Fe I is a traditional temperature 
sensor. Towards lower metallicities, higher effective tempera- 
tures, and lower gravities, the increasing UV radiation field 
and decreased rate of electron collisions may create an over- 
ionisation and over-excitation of iron with respect to the LTE 
description (iRentzsc h-Holm 1996; Mashonk ina et al.l2 01 1). The 
same effect can be found in the abundances of other neutral el- 
ements, especially those of the iron peak. Observationally, these 
non-LTE effects have been noted in, amongst others, the clas- 
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Table 4. Expected equivalent widths of the strongest He lines (at 
4471.50, 5875.65 and 5876.50 A), and iron abundances for the 
best atmosphere model for a range of effective temperatures for 
J050632. 



Table 5. Final atmosphere models. 



Teg 


Expected EW (mA) 


log e F ei° 


log e F eii fl 


(K) 


4471.50 


5875.65 


5876.50 


(dex) 


(dex) 


6250 


0.7 


0.3 


0.3 


5.69 + 0.23 


6.20 ± 0.23 


6500 


4.6 


1.1 


1.2 


5.89 ± 0.20 


6.25 ± 0.20 


6750 


10.5 


2.8 


3.0 


6.09 ±0.18 


6.28 ±0.15 


7000 


17.6 


5.6 


6.0 


6.34 ±0.17 


6.37 ±0.12 



Notes. <n) The listed errors are the line-to-line scatter. The Fel abun- 
dance is based on the measurements of 63 lines and that of Fe II on 
30. 



sical Cepheid 6 Cep (Kovtyukh & Andrievskv 1999), and some 
post-AGB stars (Takeda et al.l l2007t iGorlova et al.ll2012h . This 
means that the TefF we estimate from the excitation analyses of 
Fe I may be too high. At these temperatures this photoionisation 
scarcely alters the abundance of Fell, with differences of less 
than 0.01 dex (Mashon kina et al.1l201 ll) in the relevant temper- 
ature regime of our programme stars. The effective temperature 
deduced from a general excitation analysis of the Fe II lines is, 
however, less accurate because of the smaller number of avail- 
able features and the lower spread in excitation potentials in op- 
tical Fe II lines. 

Hydrogen line profiles do not show a strong dependence on 
metallicity of the stellar photosphere and we obtained the effec- 
tive temperature of our stars by comparing their Balmer profiles 
with the synthetic spectra of ICoelho et al.l d2005l) . When esti- 
mating the appropriate effective temperature, we focused on the 
shape of the wings, as the depth of the feature may deviate from 
the synthetic spectra. This is the case because the core of the 
line is formed in the surface layers of the photosphere where the 
standard atmosphere models are less accurate due to the occur- 
rence of surface phenomena and the boundary conditions of the 
models. We assume that the temperature and density gradient 
of our objects are similar to a model in hydrostatic equilibrium, 
despite the fact that the stars are pulsating. 

For the given effective temperature, we first determined log g 
by imposing ionisation equilibrium for iron and, when available, 
other elements with lines of diferent ions. Only lines with an 
equivalent width of less than 150 mA were used, as the equiva- 
lent width of stronger lines is no longer strongly correlated with 
the abundance of the element, because they lay on the flat part of 
the curve of growth. We tried to estimate the Fe I abundance by 
extrap olating it to an equivalen t width of mA as described in 
iKovtvukh & An drievskvl d!999l) . but the relatively large scatter 
on abundance and equivalent width of our lines rendered this im- 
possible. The average abundance was used instead. We changed 
the metallicity of the model according to the obtained Fe II abun- 
dance. Finally, we estimated the microturbulent velocity £ f in 
steps of 0.5 km/s by demanding that the abundance derived for 
the individual Fe II features does not depend on their reduced 
equivalent width. This procedure was repeated until the results 
converged towards a single model. 

The final atmosphere models for our sample can be found 
in Table [5] and their corresponding synthetic Balmer profiles in 
FigE] The final atmosphere models for both spectra of J052043 
show a slight discrepancy in preferred gravities. Because of the 
coupling of this parameter with the effective temperature, this 
indicates a small temperature difference as well. As models with 
ionisation equilibrium between Fe I and Fe II require higher val- 



Object T cS logg £ [Fe/H] 

(K) (dex) (km/s) (dex) 





6750 


0.5 


3.0 


-1.0 


_a 


5750 


0.5 


3.0 


-1.0 


_b 


5750 


0.0 


3.0 


-1.0 




5500 


0.0 


3.0 


-1.0 




4750 


2.5 


2.5 


-0.5 




A [Angstrom] 



Fig. 2. Comparison with the synthetic spectra of ICoelho et al.l 
(2005). We show for each spectrum the first four Balmer lines 
and the synthetic spectrum of the final atmosphere model (see 
also Table [5]). All synthetic spectra have been convolved with a 
Gaussian to match the resolution of the observed. 



ues of log g at higher effective temperatures, this confirms that 
J052043 became slightly cooler in December than it was in 
November. The errors on the different atmosphere parameters 
are estimated to be the grid steps and therefore no interpolation 
of models and/or new models were computed. 



4.2. Spectrum synthesis 

The abundances of some elements could not be determined on 
the basis of isolated lines. An estimate of these abundances was 
nevertheless obtained through spectrum synthesis. To match the 
synthetic line profiles with the observed ones, we need a good 
description of the line-broadening which is the combined ef- 
fect of instrumental, macroturbulent, and rotational broadening. 
Because the latter is thought to be small in F-G supergiants, 
and its effect on the spectra is difficult to distinguish from the 
macroturbulent broadening, we neglect it in our calculations. 
The macroturbulent broadening as the only unknown broaden- 
ing parameter in the line profile fitting procedure. 
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0.8 




7.5 km/s 
10.5 km/s 



W, = 42.6 mA 
log( E ) (Fe II) = 




W, = 77.9 mA 
log(£) (Fe II) = 6.27 



5413.5 5414.0 5414.5 5425.0 5425.5 

X [Angstrom] 

Fig. 3. Illustration of the necessity of different values of the 
macroturbulence at different reduced equivalent widths. We dis- 
play two Fe II lines from the spectrum of J050632 from the same 
spectral region and with a similar abundance derived from their 
equivalent width. Their macroturbulent broadening factor dif- 
fers, however, by 3 km/s. 



4.2.1. The macroturbulent broadening % m 



Based on the results iRevniers et al.l d2004l) obtained for the C- 
rich, Galactic post-AGB star IRAS 06530-0213, we expect that 
a single value for the macroturbulent broadening will not be 
sufficient to fit all lines in a spectrum equally well. This is il- 
lustrated in Fig. [3] were we show two Fell lines in the spec- 
trum of J050632 with similar abundances, as derived from their 
equivalent width, and from the same spectral region. The two 
lines clearly require a different macroturbulent broadening fac- 
tor. The reason for this dependence is not clear, but it is sug- 
gested that this is an optical depth effect where the stronger 
lines are on average being formed at lower geometrical depth 
(IRevniers et al. 2004). As similar trends have been found in 
HD 172481 ([Rgyniers & Van Winckel 1200 ll) and all other ob- 
jects in our sample, this behaviour is certainly not exceptional. 
Most authors probably do not notice this effect as usually only 
a few lines or lines in the proximity of the feature at hand are 
investigated to determine the macroturbulent broadening. 

We determined the macroturbulent broadening by a x 2 min- 
imisation of the difference between a radial-tangential broaden- 
ing profile (Grav 2008) that is corrected for the instrumental and 
microscopic broadening, and the observed feature. As the equiv- 
alent width remains unchanged while the line profile is broad- 
ened, we used the derived mean abundance of each element and 
left the macroturbulent broadening as the only free parameter to 
fit the line profile. We varied the broadening velocity in steps of 
0.5 km/s. Only the best, unblended lines of different elements 
and ionisation states were used in this procedure. 

The resulting relation between the macroturbulent broaden- 
ing and the reduced equivalent width for J050632 is shown as 
an example in Fig. [4] More information on the fitting parame- 
ters of this relation in the other spectra is listed in Table |6] The 
cl ear difference f or ioni sed and neutral species that is observed 
in lRevniers et al.l (120041) . is not recovered for our spectra and the 
slopes we find are much smaller than their slope of 10.9. 
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Fig. 4. The macroturbulent broadening £„, of J050632 as a func- 
tion of reduced equivalent width logiW^/X)- The dashed lines 
are the standard deviation on the mean trend. 



Table 6. Fitting parameters of the relation between the macro- 
turbulent broadening and the reduced equivalent width. 



Object 


Slope 


Intersection 


cf 


Corr. Coeff. 


Nr 
1 'lines 


J050632 


2.2 


20.6 


2.2 


0.60 


107 


J052043_a 


5.4 


30.4 


2.3 


0.79 


21 


J052043_b 


3.5 


25.1 


1.9 


0.78 


52 


J053250 


3.3 


22.8 


1.7 


0.64 


30 


J053253 


1.0 


8.9 


2.2 


0.26 


58 



Notes. (n) Standard deviation with respect to the fitted trend. 

4.2.2. Blending 

All the single lines used to obtain specific abundances were 
checked individually to avoid the possibility that unrecognised 
blends could cause an overestimate of the abundance. For all 
used and presumably single lines, a list of all potential neigh- 
bouring lines was obtained from the VALD (Kupka et al. 1999) 
database. Their expected equivalent widths were computed 
based on the mean abundances found for every chemical ele- 
ment. 

For each used blend, a list of all neighbouring lines in an 
interval of 1 A was extracted from VALD. We calculated the 
appropriate macroturbulent broadening with the relations from 
Sect. 14.2.11 for the total equivalent width of the line, as the el- 
ement under consideration is the dominant contributor of the 
blend. This induces a possible error, as this relation is calculated 
for unblended lines, but it is unclear how it should be adapted 
in the case of blended lines. The synth driver of MOOG and all 
known abundances were used to perform the spectrum synthe- 
sis. The final abundance was obtained by a x 1 minimisation of 
the difference between the observed spectrum and the model in 
the spectral range used to calculate the equivalent width of the 
line. The abundances of S, Na, Zn, Ba, and Sm were based on 
blended lines. 

4.3. Elements beyond the Ba-peak 

Revn iers & Van Winckell (120031) detected gadolinium, ytter- 
bium, lutetium, and maybe also tungsten in the spectra of three 
Galactic, s-process enriched post-AGB stars. Inspired by their 
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findings, we performed a specific systematic search for lines of 
elements beyond the Ba-peak. We focused on the elements from 
gadolinium (Gd, Z = 64) up to lead (Pb, Z = 82), as the im- 
portant features of elements with lower atomic numbers were 
already included in our original list (see Sect. I4.lt . 

For each element, the atomic data of al l features in the op - 
tical spectrum were acquired from VALD (Kupkaet al. 1999). 
The strongest expected lines were determined by imposing an 
ad hoc high abundance and calculating the corresponding equiv- 
alent widths for the appropriate atmosphere model. We then or- 
dered the lines according to expected relative strength and sys- 
tematically checked the obtained spectra for the presence of the 
spectral lines. 

We discarding all lines for which (1) the nearest line in the 
spectrum was not centred on the given wavelength within 1 km/s, 
(2) the signal to noise in the spectral region was too low to rec- 
ognize a genuine spectral line, or (3) the line was located in a 
crowded region where the error on the position of the continuum 
appeared too large. For the remaining features, a list of lines of 
elements with already known abundance in the interval, was ex- 
tracted from VALD and the equivalent width of these features 
determined with MOOG. Only lines of which about ninety per- 
cent of the observed equivalent width could not be explained 
by the elements of which the abundance was already known at 
that stage, were kept in the sample. After this analyses, we com- 
pared the abundances derived from all lines of the specific ele- 
ment. Features with too high an equivalent width were discarded 
as blends with unrecognised contributors. A final quality check 
was performed by recalculating the expected equivalent width 
of the strongest lines of the element at hand with this new final 
abundance, and comparing this to what is seen in the spectra. As 
the atomic data of all lines listed in VALD are not equally reli- 
able, we only executed the latter test for elements of which the 
abundance is based on at most three lines. 

The strongest features of all heavy s-process elements are 
concentrated in the low S/N region from 3280 to 4519 A and 
we were plagued by the low S/N ratio of our spectra in the blue. 
Nonetheless, we were able to determine the abundances for a 
range of heavy s-process elements. Lines of Gd(64), Tb(65), 
Dy(66), Er(68), Tm(69), Yb(70), Hf(72), Ta(73), W(74), Hg(80) 
were positively itentified in many of the objects, and the results 
are given in the abundance Tables 171 [8ll9l and [TOl 

4.4. Hyperfine structure 

We investigated the influence of hyperfine splitting (hfs = the 
splitting of spectral lines caused by nucleon-electron spin inter- 
action) on the derived abundances of all odd-Z elements from 
La (Z = 57) to Lu (Z = 71). Because hfs causes extra de- 
saturation, it will only affect the abundance determination when 
this is based on strong lines that are situated on the flat part of 
the curve of growth. 

To make sure our results are independent of the method 
used for abundance determination, we first recalculated the 
abundances of the elements by use of spectrum synthesis (see 
Sect. I4.2.2l i. treating each feature as a single line. This process 
was then repeated taking hfs into account, with the final goal of 
comparing the acquired abundances. For each spectral line, we 
obtained the energy levels corresponding to its wavelength from 
the National Institute of Standards and Technolog^Q (NIST). 
Appr opriate value s of th e hfs constants A and B were obtained 
from lLawler et all (12001 al) for La II (Z = 57), from lSneden et all 



(120091) for PrI I (Z = 59), from lLawler et all (l2001bl) for EuII 
(Z = 63), fromlLawler et al. (2001c}forTbII(Z = 65), and from 
Sneden et al. (2003) for Lull (Z = 71). We ignored Pm (Z = 61, 
no stable isotope), Ho (Z = 67, no detections) and Tm(Z = 69, 
no A,B constants found) in our calculations. Because the electric 
quadrupole interaction (B constant) has a much smaller effect on 
the line component pattern than the magnetic dipole interaction 
(A constant), it is often neglected and we set its value to zero if 
it could not be found in the literature. The relative strengths of 
the diffe rent hfs com ponents were calculated with the equations 
given in Condon & Shortley] d!935l) . 

A real, but small effect of at most A a b = log 6m s - log e no hfs = 
-0.10 dex is noticed in some features, but falls within the orig- 
inal error margin on the abundances, and is dwarfed by other 
uncertainties like e.g., the exact value loggf, the placement of 
the continuum or undetected blends. We conclude that the effect 
of hfs on the derived abundance is very small. 

5. Abundance results 

In Tables|7j[8j|9l and[10]we present the complete abundance ana- 
lysis results of the stars in our sample, the details on all lines used 
in this analysis can be found in five catalogues that are available 
at the CDS0 and contain the following information. For each ca- 
talogue, Col. 1 gives the name of the ion a certain line belongs 
to, Col. 2 contains its rest frame wavelength, Col. 3 lists its lower 
excitation potential. The logarithm of the oscillator strength can 
be found in Col. 4, Col. 5 gives the measured equivalent width, 
and the abundance deduced from this line is mentioned in Col. 6. 
We determined the errors on the [el/Fe] ab undances in 
Table s |7j |8] |9] and [TOl with the method described in iDeroo et al.l 
(2005), and the alternative procedure for estimating the error 
induced by the parameter uncertainty of the atmosphere model 
from Revniers et al. (2007). To estimate this error, we calculated 
the abundance difference for a certain element for our preferred 
model and for the most appropriate model to fit our data with 
a 250 K higher/lower temperature and 0.5 higher/lower gravity 
(log g). The total uncertainty on the [el/Fe] abundances can be 
found in Tables |7j [8] [9] and [10] It is the quadratic sum of the 
uncertainties on the mean due to line-to-line scatter (cri t i), the 
abundance ratio uncertainties induced by the model parameters 
(o"teff > °~iogg), and the uncertainty on the Fe abundance (cr Pe y. 



Ctot 




■ (cTteff) 2 + (criogg) 2 + 



0"Fe 



If less than 5 lines were available, a line-to-line scatter of 
0.2 dex was applied. All results are presented graphically in 
Figs. [5] and [6]. We now highlight the most important abundance 
characteristics. 



5.1. Metallicity 

J050632, J052043, and J053250 are strongly metal deficient with 
a [Fe/H] of - 1 .2 (see TableQI). J053253 is only mildly metal de- 
ficient with [Fe/H] = -0.5. Overall, the metallicities of the stars 
in the LMC range from [Fe/H] = -2.0 u p to [Fe/H] = -0.3 , 
with most objects having [Fe/H] ~ -0.5 dGeisleret ail 12009). 
The fact that three of the four programme stars display metal- 
licities on the low-end of the detected range, indicates that these 
are old, evolved objects with low initial masses. 



http : //www .nist . gov/pml/index . cfm 



3 via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ A/ 



Els van Aarle, Hans Van Winckel, Kenneth De Smedt, Devika Kamath, Peter R. Wood: The s-process in the LMC 
Table 8. Same as TableH but for the spectra of J052043.86-692341 .0 taken in November and December 2008. 
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For each object, the other iron peak elements follow this de- 
ficiency but the number of detected useable lines is much lower. 
The relative abundances of nickel in J050632 and of zinc in 
J052043 and J053253 are somewhat low. 

The non-LTE ionisation effect in J050632, J052043, and 
J053250 can be seen in the abundances of the other elements 
of the iron peak if results of more than one ionisation state are 
present. Both titanium and chromium show this bias, with de- 
rived abundances on the basis of neutral lines being significantly 
lower than the abundances based on ionised lines. This non- 
LTE effect is not seen is silicon and calcium where abundances 
based on differen ions are in accordance. Only in the spectrum 
of J053250 does the silicon abundances deviate from this trend 
with Si I having an abundance of 6.85, while it is 7. 12 for Si II. In 
the abundance pattern of J053253, the non-LTE effect is not wit- 
nessed, and the obtained abundance of Cr I is even higher than 
that of CrII. 



5.2. C/O-ratio 

All objects for which the determination of the carbon abundance 
was possible, are clearly enriched, with values of [C/Fe] rang- 
ing from 1.5 in J050632 to 2.5 in J053250. Surprisingly, large 
enhancement in oxygen is also observed in these objects with 
[O/Fe] around 0.8 in J052043, J050632 and J05325. The C/O 



number ratios remain relatively moderate as a consequence of 
this large oxygen enhancement, with values of 2.6 to 1.5 (see 
Table [TTIl. No conclusions on both the carbon and oxygen abun- 
dance of J053253 are possible as at the low effective temperature 
of this object, no lines of these elements are visible in the optical 
wavelength domain. 



5.3. a-elements 

The mean of the [el/Fe] ratios for the a-elements Mg, Si, 
S, Ca, and Ti, falls between 0.1 and 0.3 for all objects (see 
Table ITTb . This is slightly deficient when compared to stars of 
the same metallicities in the Galaxy, but consistent with the a- 
enhanc ement expec ted for LMC stars at their respective metal- 
licities (Pompeiaet al. 2008). Because the cr-elements are pre- 
dicted to be mainly produced in high-mass type II supernovae, 
iPompeia et alj (2008) interpret the moderate c-enhancement at 
low metallicites in the LMC as indicative of a stronger influence 
of type la supernovae. This also accords with the more extended 
star formation history in the LMC than in the Galaxy. 



5.4. s-process elements 

In general, the s-process elements observed in evolved stars are 
subdivided into three groups, depending on the number of neu- 
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Table 7. Abundance results for J050632. 10-7 14229.8. 
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Table 9. Same as TableH but for J053250.69-713925.8. 



Notes. For each element, we list the number of lines involved in the 
abundance determination, the average equivalent width of these lines, 
the deduced abundance, the standard deviation on this abundance, and 
the [el/Fe] value with its corresponding error (se e text for de t ails). T he 
solar abundances in the last column are from Asolund et alj (l2009h . If 
no solar photospheric abundance value was available, as was the case 
for Ta, Re, and Hg, we took the meteoric equivalent. 
(fl) The indication ss denotes that the abundance was deduced with the 
use of spectrum synthesis (see Sect. 14.2.21 . 

(b) These abundances should be treated with caution as they are based 
on small features when compared to the noise level, and are upper limits 
only. 



trons in their nuclei. The light s-process elements have neutron 
numbers around the magic number 50 (Y, Zr, . . . ), and the neu- 
tron numbers of the heavy s-process elements are distributed 
around the magic number 82 (Ba, La, Ce, Pr, Nd, Sm, . . . ). The 
third and final group corresponds to the double magic state of 
neutron number 126 and consists of lead only ( 208 Pb contains 
126 neutrons and 82 protons), as it is seen as the end product of 
the s-process nucleosynthesis. 

From Table QT| and Figs. and [6] it is clear that all ob- 
jects are strongly enriched in s-process elements, which con- 
firms the post third dredge-up status of our stars. The enrich- 
ment in s-process elements is the lowest in J053253, where 
[s/Fe] = 0.85. J053250 and J052043 a re s-process en hanced 
at almost the level of I RAS06530-0213 dRevniers et alj [2004), 
and IRAS 05341+0852 (IVan Winckel & Revniersl l2000h . which 
are, with respectively [s/Fe] =2.1 and [s/Fe] = 2.2, some of 
the most s-process enriched intrinsic Galactic objects known to 
date. 
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6. Luminosities and initial masses 

With the best atmosphere models as followed from the abun- 
dance analysis, we are now able to make a better estimate of the 
luminosities, and consequently initial masses of the stars in our 
sample. 



6.1. Luminosity 

In Fig. [7] we show the SEDs of the different stars, where the 
appropriate atmosphere model (see Sect. I4.lt was fitted to the 
dereddened photometry. The E(B-V) was determined by min- 
imising the difference between the scaled model and the dered- 
dened photometric data, excluding the infrared where thermal 
emission of dust contributes to the total fluxmeasurement. Errors 
in the luminosity and the reddening are dominated by the un- 
certainty in the parameters of the model atmosphere which are 
more important than the uncertainties of the photometric data 
points themselves. We therefore ran our SED fitting routine for 
all input models within the uncertainties of the spectral analyses. 
In Table fT2l we compare our luminisity estimates of the dif- 
ferent program stars, with those obtained by integrating the raw 
photometric data over the full SED (in the table given in the 
column raw SED integral). The luminosity based on integral of 
the raw data holds a good estimate of the luminosity based on 
the model atmosphere fit. This indicates that the circumstellar 
shell is spherically symmetric, and that the reddening caused by 
the interstellar medium is negligible. With these assumptions, 
the interstellar reddening can be neglected, while the extinction 
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Table 11. The C/O number ratio, metallicity, enhancement in a-elements and the s-process indexes for the stars in our sample. 
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1.75 ±0.06 


1.74 ±0.07 


J053250 


2.5 ± 0.7 


-1.22 ±0.11 


0.26 ± 0.07 


1.83 ±0.07 


1.78 ± 0.07 


J053253 




-0.54 ± 0.09 


0.18 ±0.07 


0.85 ± 0.07 


0.75 ± 0.08 


Object 


[ls/Fe] 


[hs/Fe] 


[hs/Fe] LaNd 


[hs/ls] 


[hs/ls] La-Nd 


J050632 


1.42 ±0.05 


1.06 ±0.14 


1.33 ±0.14 


-0.36 ±0.15 


-0.09 ±0.15 


J052043.a 


1.65 ±0.10 


1.95 ± 0.07 


1.84 ±0.07 


0.30 ±0.13 


0.39 ±0.13 


J052043_b 


1.58 ±0.10 


1.79 ±0.08 


1.91 ±0.09 


0.21 ±0.13 


0.34 ±0.14 


J053250 


1.53 ±0.11 


1.96 ± 0.08 


2.03 ± 0.09 


0.43 ±0.14 


0.50 ±0.14 


J053253 


0.55 ±0.10 


0.94 ±0.10 


0.95 ±0.13 


0.39 ±0.14 


0.41 ±0.16 



Table 10. Same as Table [7] but for J053253.51-695915.1. For 
this cooler star, Fe I is the dominant ionisation stage. 
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1 


47 


1.48 


0.20 
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0.36 
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Table 12. Comparison of luminosity and reddening values ob- 
tained with our different computation methods. 



caused by the circumstellar shell is accounted for by the lumi- 
nosity offlieJRexcess. 

In lvan Aarle et al.l (1201 ll) . we used the spectral type to deter- 
mine a range of possible effective temperatures and concluded 
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Fig. 5. The derived abundance patterns for J050622 (upper 
panel) and J052043, obtained with highest S/N spectrum (lower 
panel). 



simultaneously on the best reddening factor and effective tem- 
perature by fitting atmosphere models to the dereddened pho- 
tometry. We used, however, only models with metallicity -0.3, 
which is too high for three of the four stars we discuss here in 
detail. These objects are intrinsically bluer than the used models, 
which leads to an overestimation of the effective temperature, 
and consequently a too high reddening and luminosity. This lead 
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Fig. 6. The derived abundance patterns for J053250 (upper 
panel) and J053253, obtained with highest S/N spectrum (lower 
panel). 
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Fig. 7. Final SEDs of the stars in our program. We fitted the 
dereddened photometry to the rescaled atmosphere models, the 
parameters of which are based on the abundance analysis. 
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Fig. 8. The position of our stars i n the HR diagram. The dis - 
played evolutionary tracks are from lVassiliadis & Woodl (11994). 
with metallicities of Z = 0.001 (upper panel), and Z = 0.004 
(lower panel), and are labelled with their initial mass in M , The 
location of each star is shown in black in the appropriate diagram 
and in grey in the other to compare the position of the different 
data points. 



in Ivan Aarle et alj (1201 lb to an overestimate of the luminosity 
and the integral over the raw photometry is a better estimate of 
the true luminosity. In Kamath et al. (2013) we report on our 
low-resolution spectral fitting routine to obtain more accurate 
model photospheres on the basis of low -resolution spectra. 

6.2. Initial mass 

The more accurate values of effective temperature and lumi- 
nosity allow us to estimate the initial mass of the stars in our 
sample. In Fig. [8] the position of the objects in the HR di- 
agram is shown, together with the evolutionary tracks from 
Vassiliadis & Wood] d!994l) . The tracks in the upper plot have 
a metallicity of Z — 0.001 which is th e lowest value in the 
model grid of I Vassiliadis & Wood! (1 19941) but is still somewhat 
higher than what we find for J050632 (Z = 0.0007), J052043 
(Z = 0.0008), and J053250 (Z = 0.0007). The tracks in the 
lower plot correspond to Z = 0.004, which is the grid value that 
approaches the metallicity of J053253 (Z = 0.0040) the best. 
We converted the metallicities of the stars in our program to Z- 
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values by computing this parameter with its definition, under the 
assumptions that the Fe-abundance remained the same and that 
all initial abundances, except for those of hydrogen and helium, 
scale with the iron abundance. 

From Fig. [8] it follows that the initial masses of the stars in 
our sample are all smaller than about 1.5 M , under the assump- 
tion that these objects follow the single star evolutionary tracks. 
J052043 appears to be the most massive, followed by J053250 
and J050632, although the difference is small. J053253 is located 
below the evolutionary tracks, which would indicate an signifi- 
cantly lower initial mass. This is difficult to understand as ob- 
jects with masses below 0.8 M have main sequence life times 
that are larger that the Hubble time. We conclude that J053253 
is likely not evolving along the single-star evolutionary tracks 
depicted here. The evolutionary status of this object is discussed 
further in Section 8. 

7. Neutron exposure 

Four observational indexes are traditionally defined to describe 
the s-process overabundances as well as the s-process distribu- 
tion: [s/Fe], [ls/Fe], [hs/Fe], and [hs/ls]. The specific elements 
taken into account when calculating these indexes, vary from au- 
thor to author, and mainly depend on which abundances could be 
determined and the reliability of these results. To b e consistent 
with earlier publicat ions (e.g., Revniers et al.ll2004l) . we follow 
the suggestion from lBusso et al.l (11995) and define the ls-index 
as the mean of the relative abundances of Y, and Zr, and the hs- 
index as the mean of Ba, La, Nd, and Sm. Consequently, [s/Fe] 
is the mean of the six elements and [hs/ls] = [hs/Fe] - [ls/Fe]. 

The s-process indexes of our sample stars are listed in 
Table QT| Missing abundances of specific s-process elements 
need to be determined, however. For this we scaled the abun- 
dance distribution prediction of the most appropriate AGB 
model (see Section |9]i to the nearest element in atomic mass 
with a measured abundance. We then assume that the object dis- 
plays the same abundance ratio for elements nearest in atomic 
mass. The Ba abundance of J052043 _a, J052043_b, J053250, and 
J053253, and the Sm abundance of J053250 had to be estimated 
this way. The four heavy s-process elements all originate from 
the same iron seeds and we assume that their relative abundances 
only depend on their cross-sections for neutron capture. Leaving 
the missing abundances out is not a reliable alternative as due to 
the odd-even effect, the overabundances of these four elements 
are not the same. Alternative indexes of compar able quality but 
based on the abundances of La, and Nd only dRevniers et al.l 
2004) can also be obtained, and we list these for comparison. 

Observationally, it is well known that there is a strong 
correlation between th e [s/Fe] and [hs/ls] indexes (e.g., 
I Van Winckel & Revniers! l2000t IRevniers et alfl2004l) . with ob- 
jects with a higher [s/Fe] index showing also a higher [hs/ls] 
ratio. The [s/Fe] index forms an indication of the third dredge- 
up efficiency, although this parameter is also influenced by the 
mass loss history of the star and the envelope mass at dredge- 
ups. The [hs/ls] index represents the effectiveness of the neutron 
irradiation, which is related to the ratio of the number of neu- 
trons available to the iron seed nuclei. Considering that the stars 
in our sample are of similar metallicity, this parameter would be 
determined by the abundance difference between 13 C, the main 
neutron se ed, and 14 N, th e main neutron poison in the stellar in- 
terior (e.g. Herwig 2005). This difference, however, depends on 
many unconstrained parameters. The neutron production in AGB 
models is fuelled by the injection of protons into the intershell, 
which creates the 13 C pocket through 12 C(p,y) 13 N08 + ) 13 C, and 




2.5 



Fig. 9. The correlation between the total enrichment in s- 
process elements and the [hs/ls] index. The objects of this 
paper are na med. The others are the 21 //m sources from 
Van Winckel & Revniers (2000), and the two enriched objects 
from IRevniers et all (120041) . The full line gives the least-squares 
fit to all data points, the dashed line to all data points except for 
J053253, and the dotted line, which is given for comparison, to 
all older data points. 



allows the neutron producing reaction u C(a, n) 16 to take place. 
The correlation therefore indicates that with increasing dredge- 
up efficiency, also the dredge-in of protons will increase as they 
are the primary cause of a higher neutron irradiation. 

The correlation between the strength of the neutron irradi- 
ation, and the efficiency of the third dredge-up is graphically 
depicted in Fig. [9] Three of the four stars in our sample fur- 
ther confirm this assertion, but in J053253 the relatively high 
neutron irradiation only led to a mild enhancement in s-process 
elements. A simple, linear least-squares fit to all data points 
gives [hs/ls] = 0.55 x [s/Fe] - 0.78, and a correlation co- 
efficient of 0.52. If we ignore J053253, this relation becomes 
[hs/ls] = 0.92 x [s/Fe] - 1.43, with a correlation coefficient of 
0.92. This is comparable to the cor relation coefficient o f 0.96 
that corresponds to the original data (IRevniers et al.l l2004). 

The [hs/ls] index is expected to increase with decreasing 
metallicity as a consequence of the larger number of neutrons 
per iron seed nucleus. This assumes, however, that the dif- 
fusion of protons into the C-rich intershell itself is indepen- 
dent of the metallicity or parameters linked to this metallicity. 
Only a weak correlation with a broad intrinsic spread is ob- 
served between the neutron irradiation and the metal licity (e.g., 
IVan Winckel & Revniersll2000t IRevniers et al.ll2004l) . which in- 
dicates that other stellar parameters also determine the internal 
nucleosynthesis during the AGB evolution, and we are dealing 
with a large spread in 13 C pocket efficiencies in highly similar 
objects. 

With their range of metallicities from -0.54 to -1.22, our ob- 
jects offer the possibility to extend this relation to lower val- 
ues of [Fe/H]. In Fig. [TU]we a dded the current sample of post- 
AGB stars to th e objects from lVan Winckel & Revniersl fcOOO) 
and Revniers et al. (2004) in a plot of [hs/ls] versus [Fe/H]. Our 
results confirm that there is no correlation between metallicity 
and the [hs/ls] index in this range of metallicities. 
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Fig. 10. The [hs/ls] index as a function of the metallicity. The 
symbols are the same as in Fig. [9] No correlation is found be- 
tween the metallicity and the [hs/ls] index in this range of metal- 
licities. 



8. Is J053253 an extrinsic object? 

Objects that are enhanced in s-process elements, can have 
reached this status in two different ways. In the case of intrinsic 
post-AGB stars, the surface composition became s-process en- 
riched by dredge-ups of nucleosynthetic products created in the 
stellar interior. In extrinsic objects, the enrichment is an effect of 
mass transfer of s-process elements from a further evolved com- 
panion to the star at hand. The former AGB star is now typically 
a cool white dwarf. Therefore, it is possible that the enhanced 
object itself is not yet a post-AGB star. 

The strange position of J053253 in Figs. [8] and [9] suggests 
that the s-process enrichment in this object may be extrinsic: 
J053253 is much less luminous than can be expected from post- 
AGB evolutionary tracks, and it does not follow the strong corre- 
lation that is found between third dredge-up efficiency and neu- 
tron irradiation. There are, however, two arguments that contra- 
dict this interpretation. First of all, ITomkin et al.l (1989) detected 
only two barium dwarfs in their survey of about 200 F-G dwarfs, 
which is co mparable to the about on e percent of RGB stars that 
is extrinsic (MacConnel l et al.lll972l) . A similar number is hence 
expected for post-AGB stars, and it would be very surprising to 
encounter one in our sample of only four stars. 

The second argument considers the dust that is seen in the 
system of J053253. This object clearly has a dust excess, al- 
though it is rather small when compared to the other stars de- 
scribed here. These excesses are generally not seen in extrinsic 
objects, as the circumstellar dust dissipates into the interstellar 
medium (ISM) on a timescale that is short compared to the evo- 
lutionary timescales. The relative proximity of the dust in the 
system of J053253, as follows from its colour, indicates that, if 
the material was lost in a wind, it must have been shed very re- 
cently. In this scenario we may conclude that J053253 is either 
indeed an intrinsic post-AGB star, be it at very low luminosity or 
it is a post-RGB star which happens to have had a dusty mass- 
loss event. A post-RGB status conflicts with internal s-process 
enhancement. 



Another possible explanation for the deviating position of 
J053253 in Figs. |8]and|9]may be found in the shape of its SED 
(see Fig. [7). Although the dust in the system is cooler than is 
typically expected for a post-AGB object with a circumstellar 
disc, it radiates the bulk of its energy at wavelengths shorter than 
24 fim. It is hence a valid possibility that the circumstellar ma- 
terial surrounding this object resides in a disc rather than in a 
freely expanding, detached shell. It is known from the Galactic 
sample of post-AGB stars with a circumstellar disc, that the 
AGB evolution has been cut short b y a phase of strong binary 
interaction dVan Winckel et alj|2009t iGielen et al.ll201 lh . Proof 
of this theory is found in the orbits that a re detected, which 
are too small to acc ommodate an AGB star (IVan Winckel et all 
2009; iGorlova et al.ll2012l) . and the chemical comp osition of the 
dust, which consists of strongly processed silicates (Giele n et al.l 
1201 ll) . This truncation of the stellar evolution should translate 
itself in a subluminous object and should also alter the expected 
chemical composition of the photosphere: the [s/Fe] index is 
supposed to be lower, because fewer third dredge-ups can take 
place, while the [hs/ls] index will remain similar as each third 
dredge-up brings heavy and light s-process elements to the sur- 
face. 



9. Comparison with AGB nucleosynthetic models 

The current theoretical models of the internal nucleosynthesis 
and the subsequent photospheric enrichment processes of AGB 
stars are sophisticated and rely on the comple x interplay between 
mixing, nucleosynthesis, and mass loss (Bu sso etaD 119991 : 
iHerwi g 2005). In th ese models, extensive nuclear networks have 
been included (e.g . Cristall o et al.ll2009l 120111: iKarakasI 120101: 
Church et al. 2009; Goriely et al. 2008; Lugaro et al. 2012), non- 
convective mixing like mixing due to differential rotatio n and 
therm ohaline mixing have been implemented ( e.g. ^S iess et al 



2003t ISiessi l2007t IStancHffe & Jeffervl 120071: lAngelou et al 



2011), the effect of deep mix ing or ex tra mixing pro cesses ha s 
been critically evaluated (e.g. Karakas 2010; Bu sso et al1 l2010). 
and overshoot regimes have been explored in more detail to 
explain the appearance of the 13 C pocket as the main neutron 
source (e.g. IHerwig 2005). The models are not self-consistent 
and overshoot regimes must be incorporated. The proton inges- 
tion process is complex and detailed simulations at low metal- 
licities in 3 D yield clearly differe nt results as the traditional ID 
approaches (Stanclif fe et al.ll201lh . 

Tailored model computations for the individual stars as well 
as detailed comparison of the the different codes like we have 
done in De Smedt et al. (2012), fall outside the scope of this pa- 
per. The main discrepancy between the model abundances and 
the derived ones was that the predicted C/O ratio was much too 
high compared to the observed one. Moreover, for that one star 
in the SMC (J00444 1.04-732 136.4), the models fit the s-process 
distribution well, despite the fact that the absol ute abundances 
were n ot well matched. The main conclusion of iDe Smedt et al.l 
(120121) was that we found only a weak dependency of the the- 
oretical predictions on the adopted stellar evolution code. With 
the newly analysed objects we want to investigate whether these 
discrepancies between AGB nucleosynthesis predictions and the 
observed chemical patters are a systematic feature and hence 
found in more objects. 

For this general comparison, we used the nucleosyn- 
thetic AGB models from the online-database FRUITYQ (Franec 
Repository of Upgraded Isotopic Tables & Yields, Cristallo et al. 



http : //www . oa- teramo . inaf . it/fruity 
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Fig. 11. Some characteristics of the used nucleosynthetic AGB 
models. We show the C/O and 12 C/ 13 C ratio, and the [s/Fe] and 
[hs/ls] indexes as a function of the number of third dredge-ups 
the model has undergone. Each third dredge-up episode is rep- 
resented by a dot or a cross, depending on the nucleosynthetic 
model, and they are linked for clarity. The model with metal- 
licity Z = 0.001 is shown as connected black dots, while the 
models with metallicity Z = 0.003 and Z = 0.006 are repre- 
sented by grey dots and crosses respectively. The values obtained 
for the stars in our sample are overplotted as horizontal lines 
in the same colour (black,grey) as their most appropriate nucle- 
osynthetic model. J050632 is represented by a full, black line, 
J052043 by a dashed, black line, J053250 by a dotted, black line, 
and J053253 by a full, grey line. 



201 1 ). The models that are currently in the database focus on low 
mass AGB stars with initial masses between 1.5 and 3.0 M and 
metallicities Z in the range from 0.001 to 0.02. 

We used the initial mass and metallicity estimates from 
Sect. 16.21 to determine the model that approximated these val- 
ues the best for each object. All models have an initial mass 
of 1.5 M Q as this is the lowest value in the grid. As we sus- 
pect that J053253 is extrinsically enriched, we assume that the 
donor has this mass as well (see Sect. I6.21 >. We use a metallicity 
of Z = 0.001 for J050632, J052043, and J053250 as also this 
value is on the lower edge of the grid, and Z = 0.003 and 0.006 
for J053253, because the metallicity of this object falls between 
these two values. 

In Fig. QT| we show some characteristics of the models as 
a function of the number of third dredge-up episodes the star 
has undergone. The model with metallicity Z = 0.001 reaches 
the post-AGB phase after fourteen third dredge-ups, while this 
already happens after respectively ten and eight for the models 
with Z = 0.003 and 0.006. We plotted the C/O and 12 C/ 13 C ra- 
tio, as well as the [s/Fe] and [hs/ls] indexes as a function of the 
number of dredge-ups the star has undergone at that point. 

Also for these stars, we confirm that the relative abundances 
of both C and O are badly predicted by the models: C is in gen- 



eral overpredicted while the model values for O remain well be- 
low the actual determined value. This is a fundamental problem 
of the models, as especially for the C-abundances, the model 
values lie well outside the error bars on our relative abundances. 
The combination of both deviations leads to a too high predicted 
C/O ratio of the AGB models. From the first panel in Fig. fl~JJ 
it follows that the model with Z = 0.001 already estimates the 
corresponding observed C/O ratios after the occurrence of only 
one third dredge-up. The C/O ratio of the final model with this 
metallicit y overrates t he observed one by a factor ~20. As in 
De Smed t et al.l d2012l) . the high observed O-abundance, which 
also contributes to the C/O ratio, is not accounted for. 

Our spectra do not allow deduction of the 12 C/ I3 C isotopic 
fractions, because the CN violet and red systems around 3875 
and 6191 A are not detected at the S/N of our spectra and the 
temperature of our stars. Our spectral range does not cover the 
additional CN red system at 6945 A, and the wavelength re- 
gion from 1-5 /im, where several CO forests could be found. 
Determining the 12 C/ 13 C isotopic ratio is an observational chal- 
lenge in these objects but it would allow a very good test on the 
predictive power of the AGB models. 

For the objects with the highest s-process enhancement, 
J052043 and J053250, the comparison of the s-process distribu- 
tion with theoretical models matches relatively well, but the nu- 
cleosynthetic models underestimate the actual abundances of the 
light s-process elements yttrium (Y, Z = 39) and zirconium (Zr, 
Z = 40) by about 0.2 dex. The underestimation of the observed 
relative abundances of the light s-process elements leads to an 
overestimation of the actual [hs/ls] index, which forms an indi- 
cation of the neutron irradiation. In the fourth panel of Fig.[TT] it 
can be seen that the models already overestimate this index after 
the second occurrence of a third dredge-up for J052043, and af- 
ter the third for J053250. The overall enhancement in s-process 
elements, or the efficiency of the third dredge-up, on the other 
hand, is relatively well reproduced for both objects (see the third 
panel of Fig. Q~TJ after the full series of predicted 3rd dredge-up 
events. 

For the other two objects, J050632 and J053253, the compar- 
ison with the theoretical models is not as good. For J053253 the 
observed abundances are systematically too low with respect to 
the standard model (Fig. [TTJ. This is another indication that the 
s-process and the third dredge-ups are indeed inhibited in some 
way (see Sect. [8}. 

The comparison for J050632 points to a systematic differ- 
ence: the model underestimates the light s-process abundances 
while it heavily overestimates the observed abundance pattern 
for the heavy s-process elements. This results in a negative 
[hs/ls] index, which is not predicted by the models (see Fig. [TTJ. 

Although J050632, J052043 and J053250 are intrinsically 
enriched post-AGB objects with a very similar metallicity, lumi- 
nosity, and hence initial mass, they display quite different levels 
of overabundance as well as a significantly different s-process 
distribution. It therefore follows that the AGB models have to be 
individually adapted to match the observed abundance patterns. 



10. Conclusions 



In Ivan Aarle et al.l d201ll) . we classified the objects J050632, 
J052043, J053250, and J053253 (see Table [TJ as spectroscopi- 
cally confirmed post-AGB stars, because of their spectral type, 
IR colours, and estimated luminosity. The chemical analysis 
we carried out in this paper was based on high-resolution 
U VES spectra, and confirmed their post-AGB status, as the third 
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dredge-up clearly took place in all objects: they are C-rich, 
and enhanced in s-process elements. We deduced abundances of 
heavy s-process elements for all stars in the sample, and quanti- 
fied also abundances of elements well beyond the Ba-peak. The 
low S/N of the spectra in the blue spectral domain, where the 
lines of heavy s-process elements abondan, prevented very ac- 
curate determinations. 

The metallicity of all stars except J053253 is considerably 
lower than the average value that is observed for the LMC. This 
is in accordance with the possibility that these objects are old, 
which means that they should have a low initial mass. The lat- 
ter is also confirmed by their low luminosities. We estimate the 
initial mass of the three objects to be between 1 and 1.5 M . 

All sample stars except J053253 confirm the correlation be- 
tween the efficiency of the third-dredge up and the neutron expo- 
sure that is present in the Galactic sample of post-AGB objects. 
The non-existence of a correlation between metallicity and neu- 
tron irradiation is also further corroborated. 

The subluminosity of J053253, its small but warm infrared 
excess, and the deviation of the s-process abundance pattern 
from the expected trend (see Fig. |9]l are all observational indi- 
cations that this object is not evolving on a single-star evolu- 
tionary track. We argued that the dust in the system of this star 
may reside in a disc rather than in an freely expanding, detached 
shell. This would mean that this object is likely a binary. From 
Galactic binary post-AGB stars with a disc, it is known that the 
AGB evolutionary p hase in these objects is tr uncated, because of 
the detected orbits dVan Winckel et all 2009). This can explain 
the subluminosity and low [s/Fe] index that are observed for this 
star, but radial velocity time-series are badly needed. 

The 21 yum feature ha s been detected in the spectrum of 
J052043 dVolk et alJl20TTh . As the fifteen Galactic objects that 
display this feature are post-ca rbon star s which are strongly en- 
hanced in s-process elements (Hrivnaket al. 2008), this detec- 
tion shows that this is also valid for the L MC. Similar LMC 
21 /im sources are known dVolket al.ll201ll) . and it would be 
useful to obtain abundances for them as well. In the SMC, 
there is on ly one known 21 //m source: J004441. 04-732136. 4 
dVolket all 1201 lh. which was recently subjected to a detailed 
chemical study dDe Smedt et alj |2012). This object has an ini- 
tial mass of~1.3 M , a metallicity of -1.3, and turned out to 
be one of the most s-process enriched objects known to date. 
A more detailed comparison with model predictions using two 
different independent codes shows that very similar conflicts ex- 
ist between the deduced and predicted abundances as found in 
this contribution. 

From the comparison with AGB models, we conclude that 
the available AGB models overestimate the observed C/O ratios 
and fail to reproduce the variety of [hs/ls] indexes that is ob- 
served in other wise very similar star s. This seems to be a com- 
mon feature dDe Smedt et al.l 120121) . Additional near infrared 
spectra would be useful to determine the 12 C/ 13 C ratio of our 
stars, as this ratio is predicted to be very dependent on the num- 
ber of dredge-up events the star has suffered on the AGB. The Pb 
production is predicted to be very high in these metal poor ob- 
jects but an observational constraint will require additional high 
S/N blue spectra. 

This paper is the first systematic study of enriched post-AGB 
stars in the LMC. Because the range of luminosities of this sub- 
sample is rather small, more general and systematic confronta- 
tions of the observed abundances with AGB nucleosynthetic mo- 
dels will only become possible when the abundances and char- 
acteristic parameters of other objects covering a wider range of 
metallicities and luminosities have been studied. We therefore 



initia ted a large systematic low-resolution survey (Kamat h et al.l 
2013) to identify the best suitable candidates both in the LMC 
and SMC which will enable us to quantify a wide range of chem- 
ical species in post-AGB stars covering a wide luminosity range. 
Despite the small range in luminosities, the detected abun- 
dances cover a range in s-process overabundances and abun- 
dance patterns which call for very detailed, tailored models that 
will match the derived abundance patterns individually. Only 
such a systematic study can make significant progress in the 
understanding of the complex interplay between mixing, nucle- 
osynthesis and mass-loss which characterise the final evolution 
of solar-like stars, and is the ultimate goal of our research. 
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